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Divergencies appearing in perturbation expansions of interacting many- 
body systems can often be removed by expanding around a suitably chosen 
renormalized (instead of the non-interacting) Hamiltonian. We describe 
such a renormalized perturbation expansion for interacting Fermi systems, 
which treats Fermi surface shifts and superconductivity with an arbitrary 
gap function via additive counterterms. The expansion is formulated explic- 
itly for the Hubbard model to second order in the interaction. Numerical 
solutions of the self-consistency condition determining the Fermi surface 
and the gap function are calculated for the two-dimensional case. For the 
repulsive Hubbard model close to half-filling we find a superconducting 
state with d-wave symmetry, as expected. For Fermi levels close to the van 
Hove singularity a Pomeranchuk instability leads to Fermi surfaces with 
broken square lattice symmetry, whose topology can be closed or open. For 
the attractive Hubbard model the second order calculation yields s-wave 
superconductivity with a weakly momentum dependent gap, whose size is 
reduced compared to the mean-field result. 
PACS: 71.10.Fd, 71.10.-w, 74.20.Mn 



I. INTRODUCTION 

Unrenormalized perturbation expansions of interacting electron systems around the 
non-interacting part of the Hamiltonian are generally plagued by infrared divergencies. 
Some of the divergencies are simply due to to shifts of the Fermi surface, while others 
signal instabilities of the normal Fermi liquid towards qualitatively different states, such as 
superconducting or other ordered phases. This problem is often treated by self-consistent 
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resummations of Feynman diagrams, where a finite or infinite subset of skeleton diagrams, 
with the interacting propagator G on internal lines, is summed.0 Symmetry breaking can 
be built into the structure of G as an ansatz, and the size of the corresponding order 
parameter is determined self-consistently. This standard approach has been very useful 
in many cases. However, resummation schemes beyond first order (Hartree-Fock) require 
extensive numerics, since the full self-energy has to be determined self-consistently, and 
delicate low-energy structures cannot always be resolved. A more serious problem is 
the fact that self-energy and vertex corrections are not treated on equal footing in most 
feasible resummation schemes. This often leads to unphysical results. 

In this work we will describe and apply an alternative procedure, which has been for- 
mulated already long ago by Nozieres,@ and more recently been discussed in the mathe- 
matical literature as a way of carrying out well-defined perturbation expansions for weakly 
interacting Fermi systems.!'! The basic idea is to choose an improved starting point for 
the perturbation expansion, by adding a suitable counterterm to the non-interacting part 
of the Hamiltonian, and subtracting it from the interaction part. The counterterm is 
quadratic in the Fermi operators and has to be determined from a self-consistency condi- 
tion. In Sec. II we will describe how Fermi surface deformations and superconductivity can 
be treated by this method. Explicit expressions up to second order in the interaction are 
derived for the case of the Hubbard model in Sec. III. Results obtained from a numerical 
solution of the self-consistency equations in two dimensions will follow in Sec. IV. For the 
repulsive Hubbard model we have obtained superconducting solutions with d-wave sym- 
metry in agreement with widespread expectations,! and with recent renormalization group 
calculations which conclusively established d-wave superconductivity at weak coupling.! 
In addition, for Fermi levels close to the van Hove singularity, deformations which break 
the square lattice symmetry occur. This confirms the recently proposed possibility of 
symmetry-breaking Fermi surface deformations (" Pomeranchuk instabilities"). EH] 



II. RENORMALIZED PERTURBATION EXPANSION 

We consider a system of interacting spin-^ fermions with a Hamiltonian H = Hq + Hj, 
where the non-interacting part 

H = ^2&n ka (1) 

k,cr 

with £k = — A* contains the kinetic energy and the chemical potential, while Hi is a 
fermion-fermion interaction term. We are particularly interested in lattice systems, for 
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which the dispersion relation ek is not isotropic. We consider only ground state properties, 
that is the temperature is zero throughout the whole article. 

The bare propagator in a standard many-body perturbation expansion^ around H 
is given by 



where ui is the Matsubara frequency and k = (a>,k). This progagator diverges for uj — > 
and k — > kp, for any Fermi momentum kp, since £k F = 0. As a consequence, many 
Feynman diagrams diverge. A well-known singularity is the (usually) logarithmic diver- 
gency of the 1-loop particle-particle contribution to the two-particle vertex in the Cooper 
channel, which leads to a (log)™ divergency of the n-loop particle-particle ladder diagram. 
This signals a possible Cooper instability towards superconductivity. Much stronger di- 
vergencies occur in diagrams with multiple self-energy insertions on the same internal 
propagator line, leading to non-integrable powers of Go (&).§§ These singularities are due 
to Fermi surface shifts generated by the interaction term in the Hamiltonian. 

The divergency problems and the superconducting instability can be treated by split- 
ting the Hamiltonian in a different way, namely as§ 



where H$ = Hq + 5 Ho and Hj = Hj — 5Hq, and expanding around Hq. The counterterm 
SH must be quadratic in the creation and annihilation operators to allow for a straight- 
forward perturbation expansion based on Wick's theorem. It is possible to chose SHq such 
that Hi does not shift the Fermi surface corresponding to H any more, and divergen- 
cies due to self-energy insertions are removed. In the superconducting state spontaneous 
symmetry breaking can be included already in 5Hq, with an order parameter Ak whose 
value on the Fermi surface is not shifted by Hj. We will now describe this procedure in 
more detail. 



G (k) 
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(2) 



iu - £k ' 



H = H + Hj , 



(3) 



A. Normal state 



A counterterm 5H = a <5£k ^ko- leads to a renormalized dispersion relation £k 
£k + <5£k in the unperturbed part of the Hamiltonian, 




(4) 
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and correspondingly to a new bare propagator 



G (k) 



1 



(5) 



The Fermi surface T associated with H is given by the momenta satisfying the equa- 
tion £k = 0. The Fermi surface of the interacting system is given by the solutions of the 
equation G _1 (0,k) = 0. This surface coincides with the unperturbed one, corresponding 
to Hq, if the renormalized self-energy E = Gq 1 — G^ 1 vanishes on JF, that is if 



This imposes a self-consistency condition on the counterterms which can be solved iter- 
atively. For isotropic systems the shift of £k can be chosen as a momentum independent 
constant, which may be interpreted as a shift of the chemical potential. For anisotropic 
systems, however, one generally has to adjust the whole shape of the Fermi surface. That 
this procedure really works at each order of the perturbation expansion has been shown 
rigorously for a large class of systems.0 

The shift function <5£k is uniquely determined by the self-consistency condition only 
on the (interacting) Fermi surface T . For momenta away from the Fermi surface, 5£k can 
be chosen to be any sufficiently smooth function of k which does not lead to artificial 
additional zeros of £k- 

The perturbation expansion of the renormalized self-energy E involves two types of 
vertices: the usual two-particle vertex given by the interaction Hi and a one-particle 
vertex due to the counterterm —SH in Hj. In Fig. 1 we show the Feynman diagrams 
contributing to E up to second order in the interaction. 

The above-mentioned divergencies of Feynman diagrams with self-energy insertions 
on internal propagator lines are removed in the renormalized expansion around H Q , since 
in products GqT,Go . . . EGo only one simple pole at k — (0, kp) survives, all other poles 
being cancelled by the corresponding zeros of the self-energy E. 



To treat superconducting states we also add counterterms containing Cooper pair 
creation and annihilation operators, in addition to a shift of £k- We consider only spin 
singlet pairing, but triplet pairing could be dealt with analogously. We thus expand 
around a BCS mean-field Hamiltonian 



E(0,k) = for kG J. 



(6) 



B. Superconducting state 




A k a_ ki a kT 



+ A k a kT a_ ki 



(7) 



k,cr k 
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where Ak is the gap-function, which has to be determined self-consistently. In terms of 
Nambu operators 

*k = ( J and *+ = (a kT , a_ ki ) (8) 

one can rewrite Hq in more compact form as 

Ho = J2 & *k *k ( A k ^ - A k "2) , (9) 

k k 

where o"i,o"2,o"3 are the Pauli matrices, and A k (A k ) is the real (imaginary) part of Ak. 
The two expressions (7) and (9) for H differ by the constant (c-number) Xlk^k; which 
must be taken into account only when absolute energies are computed. The bare Nambu 
matrix propagator Go = — {^^^g following from Hq is given by 

G.- 1 W=f i 7/" . V (10) 

y A* iuj + ^ k / 

Extending the self-consistency condition for the normal state, we now require that the 
matrix self-energy S = Gq 1 — G _1 vanishes on the Fermi surface (defined by £k = 0), 
that is 

S(0,k) = for keJ. (11) 

Thus, for uo = and k on the Fermi surface, neither the diagonal nor the off-diagonal 
elements of G ~ 1 (A;) are shifted by the interaction term Hj. The Feynman diagrams in 
Fig. 1 apply also to the superconducting case, if lines are interpreted as Nambu matrix 
propagators. 

The above renormalized perturbation theory is reminiscent of the perturbation the- 
ory for symmetry broken phases formulated by Georges and Yedidia,0 where an order- 
parameter-dependent free energy function is constructed by adding Onsager reaction 
terms to the mean field contributions, and the actual order parameter is determined 
by minimizing this free energy. 

III. APPLICATION TO THE HUBBARD MODEL 

In this section we derive explicit expressions for the self-energy and the counterterms 
for the ground state (T = 0) of the Hubbard model, up to second order in the renormalized 
perturbation expansion. 
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The one-band Hubbard modelEl 

ij CT j 

describes lattice electrons with a hopping amplitude £y and a local interaction £/. Here 
cj CT and Ci a are creation and annihilation operators for electrons with spin projection a on 
a lattice site i, and n- ia = c i(T c\ a . Note that we have included the term fiN with the total 
particle number operator N in our definition of H . The non-interacting part of H can be 
written in momentum space as Hq = £k where £k = e k — A 4 an d is the Fourier 
transform of £y. 

Our numerical results will be given for the Hubbard model on a square lattice with 
a hopping amplitude — t between nearest neighbors and a much smaller amplitude —t' 
between next-nearest neighbors. The corresponding dispersion relation is 

ek = — 2t(cos k x + cos k y ) — At' cos k x cos k y . (13) 

We now derive expressions for the self-energy and the resulting self-consistency equations 
up to second order in U. 



A. Normal state 

1. First order 

To first order in U the self-energy is obtained as 

£«(» = [/ / Go(k')e^' 0+ - 6&, (14) 
Jk' 

where J k is a short-hand notation for the frequency and momentum integral, including the 
usual factors of (2tt)^ 1 for each integration variable. The first term results from diagram 
(la) in Fig. 1, the second from diagram (lb). Note that the tadpole diagram (la) yields a 
k-independent contribution, since the Hubbard interaction is local. The self-consistency 
condition (6) for X^ 1 ) yields, after carrying out the u/- integration, 

8& = uJ 0- d 6(/i - e k , - 8b) , (15) 

to be satisfied (at least) for k G T . Since the right hand side of this condition is a 
constant, it is natural to define <5£k by this constant for all k. Using Luttinger's theorem 
one can identify the above momentum integral with the particle density per spin, such 
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that <5£k = Un/2, where n is the total density. The self-consistency condition thus yields 
the n(n) relation of the interacting system. Since the counterterm can be chosen k- 
independently at first order, it may be interpreted as a shift of the chemical potential. 



2. Second order 



The diagrams (2b) and (2c) from Fig. 1 obviously cancel each other to the extent that 
the first order diagrams (la) and (lb) cancel. Writing 5£k — + with 5^ given 
by the constant on the right hand side of Eq. (15), such that 5^ is of order U 2 for all 
k, the sum of contributions from (2b) and (2c) is of order U 3 and can thus be ignored 
at second order. Hence, only diagram (2a) contributes to the second order self-energy. 
Using the Feynman rulesS one obtains 

±^(k) = U 2 f Uo(q) G (k - q) , (16) 

J q 

where ITo^) = — j k , Go(k') Go(k' + q). Adding first and second order terms, one arrives 
at the second order self-consistency condition 

&& = uj + S (2) (0,k) . (17) 

The counterterm 5£k has to be chosen such that the above equation is satisfied for all 
k G J, that is for all k satisfying £k = 0. Since S^ 2 ^(0,k) is momentum dependent, 
cannot be chosen constant any more. As a consequence, the Fermi surface of the 
interacting system will be deformed by interactions, even if the volume of the Fermi sea 
is kept fixed. Luttinger's theorem can be used to determine the density from the volume 
of the Fermi sea as n = 2 J -d^ ©(— fk)- 



B. Superconducting state 

For the matrix elements of the Nambu propagator we use the standard notation 

1 F*{k) -G{-k) 1 

and the analogous expression for G (k). The matrix elements of the self-energy are 
denoted by 

E W =|'^> ?W |. (19) 
1 S*{k) -E(-Jfe) 1 
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1. First order 



In the presence of an off-diagonal counterterm A k , the diagonal part of S is still given 
by Eq. (14) to first order, where Go(k) now depends on the gap function: 

G (k) = . (20) 

The first order self-consistency relation (15) thus generalizes to 

/r\ d V 1 
(2^2^-^')' (21) 



with Ek = y £ k + l^kl- Note that the above integral is the BCS formula for the average 
particle density per spin. 

The off-diagonal matrix element of S is obtained from diagrams (la) and (lb) in Fig. 
1 as 

SW(k) = -U [ Fo(k') + A k (22) 
Jk' 

to first order in U, with 



u 2 + ii + |A k | 2 



The off-diagonal part of the self-consistency condition (11) follows as 

TT f d d k' A k / . 

^ = - u Jwr^- (24) 

Extended as a condition for all k (and not just on J 7 ) this is nothing but the gap equation 
for the Hubbard model as obtained by standard BCS theory. The self-consistency relation 
requires that A k be constant on the Fermi surface, such that one naturally chooses a 
constant A k = A as an ansatz for all k. A non-trivial solution A ^ of this gap equation 
can obviously be obtained only for the attractive Hubbard model (U < 0). 

2. Second order 

The diagrams (2b) and (2c) cancel each other for the same reason as in the normal 
state. The contribution from diagram (2a) to the diagonal part of the self-energy is still 
given by formula (16), with G from Eq. (20) and 

note) = - / \Go(k') G (k' + q)+ F (k') F*(k' + qj\ . (25) 
Jk' 1 J 



8 



The second order contribution to the off-diagonal matrix element of S is 



SW(k) = U 2 [ U (q) F (k - q) 



(26) 



Jq 



The self-consistency relations read: 



i£k = E/ 



/S^^ i -«-' /Bk ') +i:(2,(o - k) - 



(27) 





In Appendix A we present more explicit expressions for E( 2 )(0, k) and S^{0, k), obtained 
by carrying out the frequency integrals. 



C. Numerical solution 



The self-consistency conditions are non-linear equations for the counterterms <5£k and, 
in the superconducting state, Ak- The Fermi surface of the interacting system, J 7 , on 
which the self-consistency conditions must be satisfied, is not known a priori. The equa- 
tions involve one momentum integral at first order, and two momentum integrals at 
second order. Such a non-linear system can only be solved iteratively. In this subsection 
we describe some details of our algorithm. 

Since the counterterms are determined by the self-consistency conditions only on the 
Fermi surface, their momentum dependence away from T can be parametrized in many 
ways. We have chosen <5£k and Ak as constant along the straight lines connecting the 
square shaped line defined by the condition \k^\ + \k y \ = n with the points (0, 0) and (it, ir) 
of the Brillouin zone, respectively (see Fig. 2). For a numerical solution the remaining 
tangential momentum dependence is discretized by up to 256 points. 

The iteration procedure starts with a tentative choice of counterterms. To be able 
to reach a symmetry broken solution one usually has to offer at least a small symmetry 
breaking counterterm in the beginning.0 In each iteration step new counterterms are 
determined via Eq. (17) in the normal state, and by Eqs. (27) and (28) for the supercon- 
ducting state. The right hand side of these equations is evaluated using the counterterms 
obtained in the previous step, and k is chosen on the Fermi surface defined by the previous 
<5£k- The momentum integrals are carried out using a Monte-Carlo routine. The iteration 
is continued until convergence is achieved, that is until the counterterms remain invariant 
within numerical accuracy from step to step. In all cases studied different choices of initial 
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counterterms lead to the same unique solution. The symmetry breaking terms are much 
larger than the stochastic noise from the Monte-Carlo routine in all results shown. 

The density is kept fixed by adjusting the chemical potential during the iteration 
procedure. To avoid a higher numerical effort we have computed the density from the 
Fermi surface volume in the normal state (justified by Luttinger's theorem), and from 
the BCS formula for the density in superconducting solutions. The latter reduces to 
the Fermi surface volume in the normal state limit, such that the potential error of this 
approximation is very small as long as the gap is small. 

IV. RESULTS 

We now discuss the most interesting results obtained within the renormalized pertur- 
bation theory described above, focussing mainly on the repulsive Hubbard model (U > 0), 
for which we have found superconducting solutions with d-wave symmetry, as well as 
symmetry-breaking Fermi surface deformations. 

A. Repulsive Hubbard model 

The following results for the repulsive Hubbard model have been computed for the 
parameters t' = — O.lbt and U = 3t. The interaction is thus in the weak to intermedi- 
ate coupling regime. For too small [/-values it becomes very hard to resolve the small 
superconducting gap in the numerical solution. 

We have solved the self-consistency equations for various densities ranging from n = 
0.88 to n — 0.90, for which the Fermi surfaces are quite close to the saddle points of 
the bare dispersion relation ek, located at (tt, 0) and (0, 7r). In all cases the normal state 
is unstable towards superconductivity. The gap function in the superconducting state 
obtained from the self-consistency equations has d x 2_ y 2-wave shape, with slight deviations 
from perfect d-wave symmetry in cases where the Fermi surface breaks the symmetry of 
the square lattice. This is in agreement with widespread expectations for the Hubbard 
model,! and in particular with recent renormalization group arguments and calculations.! 
In Fig. 3 we show the gap functions obtained at the densities n = 0.88 and n = 0.9, 
respectively. We note that the size of the gap is roughly one order of magnitude smaller 
than the critical cutoff scale A c at which Cooper pair susceptibilities diverge in 1-loop 
renormalization group calculations for comparable model parameters.© There are various 
possible reasons for this quantitative discrepancy. First, and probably most importantly, 
the enhancement of effective interactions due to fluctuations, especially antiferromagnetic 
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spin fluctuations, is captured much better by a renormalization group calculation. Second, 
the approximate Fermi surface projection of vertices driving the renormalization group 
flow can lead to an overestimation of effective interactions and hence of critical energy 
scales. Furthermore, a renormalization group calculation within the symmetry broken 
phase could yield a gap that is somewhat smaller than A c . 

While superconductivity is the only possible instability of the normal Fermi liquid state 
in the weak coupling limit (except for the case of perfect nesting at half-filling), at higher 
U one should also consider the possibility of other, in particular magnetic, instabilities. 
This could be done within renormalized perturbation theory by allowing for counterterms 
introducing magnetic or charge order. 

The Fermi surface is always deformed by interactions. The shifts generated by the 
momentum dependence of the counterterm 5^ are not very large. They are more pro- 
nounced near the saddle points of ek, where small energy shifts lead to relatively large 
shifts in k-space. However, the results presented in Fig. 4 show that the Fermi surface 
of the interacting system can nevertheless differ strikingly from the bare one. For the 
densities n = 0.88 — 0.889 the Fermi surface of the interacting system obviously breaks 
the point group symmetry of the square lattice. For n = 0.88 and n = 0.888 even the 
topology of the Fermi surface is changed by interactions. The deformed surface has open 
topology in these cases, instead of being closed around the points (0, 0) or (tt, it) in the 
Brillouin zone. Note that the symmetry-broken Fermi surfaces shown here correspond to 
stable solutions of the self-consistency equations for the counterterms, while symmetric 
solutions are unstable. 

More details about the Fermi surface shifts can be extracted from a plot of the second 
order counterterms, shown in Fig. 5. The actual shifts are determined by these terms plus 
a constant due to the first order counterterm and a shift of the chemical potential. At fixed 
density the interaction shifts the Fermi surface outwards at points where S(0, kp) has an 
absolute minimum, and inwards at points corresponding to absolute maxima. Interactions 
thus reduce the curvature of the Fermi surface near the diagonals in the Brillouin zone. 
Fig. 5 reveals that the Fermi surface deformation is slightly asymmetric also for n = 0.9, 
but the symmetry breaking is too small to be seen in Fig. 4. 

If the Fermi surface breaks the square lattice symmetry, the gap function cannot 
have pure d-wave symmetry any more. See, for example, the gap function at density 
n = 0.88 in Fig. 3. The deviation from perfect d-wave form is however quite small, since 
the symmetry breaking Fermi surface deformation is small. 

Interaction-induced Fermi surface deformations which break the symmetry of the 
square lattice have already been discussed earlier in the literature. Yamase and Kohnoi 
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have obtained symmetry-broken Fermi surfaces within a slave boson mean-field theory 
for the t-J model. The effective interactions obtained from 1-loop renormalization group 
flows for the Hubbard model also favor symmetry-breaking Pomeranchuk instabilities of 
the Fermi surface, if the latter is close to the van Hove points.@ A systematic stability anal- 
ysis of the Hubbard model using Wegner's Hamiltonian flow equation method confirmed 
that symmetry breaking Fermi surface deformations are among the strongest instabilities.i 
It remained an open question, however, whether such Fermi surface instabilities would be 
cut off by the superconducting gap. We have observed within our renormalized perturba- 
tion theory that symmetry breaking Fermi surface deformations occur indeed more easily, 
if the system is forced to stay in a normal state, by setting A^ = 0. Whether a symme- 
try broken Fermi surface and superconductivity coexist can be seen only by performing 
a calculation within the symmetry-broken state. This has not yet been done using the 
renormalization group or flow equation methods. 

From a pure symmetry-group point of view the symmetry breaking generated by the 
Pomeranchuk instability is equivalent to that in "nematic" electron liquids, first discussed 
by Kivelson et alH. These authors considered doped Mott insulators, that is strongly 
interacting systems. A general theory of orientational symmetry-breaking in fully isotropic 
(not lattice) two- and three-dimensional Fermi liquids has been reported by Oganesyan et 
alJll Superconducting nematic states, in which discrete orientational symmetry breaking 
develops in addition to d-wave superconductivity, have been considered recently by Vojta 
et al.0 Motivated by experimental properties of single-particle excitations in cuprate 
superconductors they performed a general classification and field-theoretic analysis of 
various phases with an additional order parameter on top of d x 2_ y 2-pairing. 



B. Attractive Hubbard model 



For the attractive Hubbard model (U < 0) the renormalized perturbation expan- 
sion yields s-wave superconductivity already at first order, which is equivalent to BCS 
mean-field theory.^ At this level the gap function is constant in k-space. Extending the 
calculation to second order, a weak momentum dependence of Ak is generated, as seen in 
Fig. 6 for the parameters U = —2t, t' = — 0.15t and n = 0.9. More importantly, the over- 
all size of the gap is strongly reduced by fluctuations included in the second order terms. 
The average gap in Fig. 6 is only one third of the corresponding mean-field gap. It has 
been pointed out previously that fluctuations not contained in mean-field theory reduce 
the size of magnetic and other order parameters even in the weak coupling limit .00 
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V. CONCLUSION 



In summary, we have formulated a renormalized perturbation theory for interacting 
Fermi systems, which treats Fermi surface deformations and superconductivity via addi- 
tive counterterms. This method is very convenient for studying the role of fluctuations 
for spontaneous symmetry breaking in a controlled weak-coupling expansion. A concrete 
application of the expansion carried out to second order yields several non-trivial results 
for the two-dimensional Hubbard model. In particular, for the repulsive model we have 
obtained the gap function of the expected d-wave superconducting state and, for Fermi 
levels close to the van Hove energy, an interacting Fermi surface with broken lattice sym- 
metry, and in some cases even open topology. The symmetry-breaking pattern of the 
states with symmetry-broken Fermi surfaces is equivalent to that of "nematic" electron 
liquids discussed already earlier from a different point of view.00 

The present work can be extended in several interesting directions. After fixing the 
counterterms one can compute the full momentum and energy dependence of the self- 
energy, and hence the spectral function for single-particle excitations. At second order the 
combined effects of symmetry breaking and quasi-particle decay are captured. Allowing 
for other symmetry-breaking counterterms, for example spin density waves, one can study 
the competition of magnetic, charge, and superconducting instabilities, as well as their 
possible coexistence. Finally, the formalism can be extended to finite temperature. In 
that case the singularities of the bare propagator are cut off by the smallest Matsubara 
frequency, but Fermi surface shifts and symmetry breaking can still be conveniently taken 
into account by counterterms. 
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APPENDIX A: FREQUENCY INTEGRALS 

The Matsubara frequency integrals in the second order self-energy contributions can 
be carried out analytically by using the residue theorem. We only present the results for 
the superconducting case; the normal state results can be recovered by setting A k = in 
the following expressions. 

The frequency integrals relevant for the evaluation of n defined by Eq. (25) are 
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/ 



r (h\ r (h _L n\ - — + ^ k+q ^Ck+q - ^k^k+q 

— Go(Jb) G (k + q) - 2EkEk+q go2 + [ £ k + W 



and 



/ 



^0 Ck-^k+q - -^k£k+q /..n 

2£ k £ k+q g 2 + [£ k + £ k+q ] 2 1 ' 



t^+^ ^r^w - <A2) 



The imaginary part of Ilo does not contribute to S(0, k). Carrying out the go-integral in 
Eqs. (16) and (26) yields: 



t^(0,k) = -U 2 / / f k _ q C(k,k',q) (A3) 



q Jk' 



where 



S^(0,k) = U 2 [ [ A k _ q C(k,k',q), (A4) 



C(k, k', q) = k+q ^k+q k + 

4-Ek-q-Kk'-^k'+q |_-^k-q + -^k' + -^k'+qj 
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FIGURES 
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2a 2b 2c 

FIG. 1. The Feynman diagrams contributing to the renormalized self-energy S at first and 
second order perturbation theory; the two-particle vertices represent the antisymmetrized in- 
teraction, one-particle vertices the counterterm, and lines the renormalized bare propagator 
G . 
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FIG. 2. Parametrization of counterterms in the first quadrant of the Brillouin zone; the 
counterterms are constant along the straight lines connecting the line from (it, 0) to (0, it) with 
the points (0,0) and (tt,tt), respectively; the dashed line illustrates a typical Fermi surface. 
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FIG. 3. Gap function for n = 0.88 (larger amplitude) and n = 0.9 (smaller amplitude) as a 
function of the angle with respect to the A^-axis. 








FIG. 4. Fermi surfaces of the interacting system for different densities n. 
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FIG. 5. Second-order counterterms S^(0, kp) as a function of the angle with respect to the 
fer-axis, for the densities n = 0.888,0.889,0.9 (from bottom to top). 
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FIG. 6. Gap function for n = 0.9 as a function of the angle with respect to the /c^-axis for 
the attractive (U = —2t) Hubbard model. 
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